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The generic square-gradient density-functional model with triple-parabolic free energy is used 
to study the stability of a cavity introduced into the stretched liquid. The various properties of 
the critical cavity, which is the largest stable cavity within the liquid, are compared with those of 
the critical bubble of the homogeneous bubble nucleation. It is found that the size of the critical 
cavity is always smaller than that of the critical bubble, while the work of formation of the former 
is always higher than the latter in accordance with the conjectures made by Punnathanam and 
Corti [J. Chem. Phys. 119, 10224 (2003)] deduced from the Lennard- Jones fluids. Therefore their 
conjectures about the critical cavity size and the work of formation would be more general and valid 
even for other types of liquid such as metallic liquid or amorphous. However, the scaling relations 
they found for the critical cavity in the Lennard- Jones fluid are marginally satisfied only near the 
spinodal. 
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I. INTRODUCTION 

The bubble plays a crucial role in many natural pro- 
cesses as well as in industrial practice or even in some in- 
dustrial hazards^^. Therefore, one of the basic mecha- 
nisms of the formation of the bubble called homogeneous 
bubble nucleation (cavitations) has attracted continuous 
attention for many years from a fundamental point of 
view as well as from technological interests^£&£. The 
homogeneous nucleation occurs within a bulk fluid that 
is under-saturated and is held at a pressure lower than 
its coexisting vapor pressure at the given temperature. 
Then the fluid is said to be stretched. Another mech- 
anism of the formation is the heterogeneous nucleation, 
where the bubble formation is induced by impurities in 
the liquid or by the wall of the container of the liquid. 

A tiny embryo of bubble called the nucleus forms 
within the stretched fluid from thermal fluctuation. The 
work of formation W of such a tiny bubble is positive, 
and has a maximum W* as the function of the size or 
radius of the embryo. The nucleus with this critical size 
corresponding to the maximum W* is called the critical 
nucleus. Therefore, a bubble that is smaller than this 
critical size is unstable and will shrink and disappear im- 
mediately. As soon as the size of the embryo becomes 
larger than this critical size, it starts to increase its size 
as the larger bubble becomes more favorable energeti- 
cally. Thus the formation of bubbles is the activation 
process. The central quantities of the nucleation phe- 
nomena is the nucleation rate J, which is the number 
of critical nuclei formed per unit time per unit volume. 
Usually it is written in Arrhenius form 

^exp(-i^), (1) 

where A is a kinetic pre-exponential factor which is be- 
lieved to be weakly dependent on temperature T, and 



/cb is the Boltzmann's constant. Therefore, the temper- 
ature dependence of the nucleation rate is controlled by 
the work of formation W*. If any impurities or walls 
exist this will assist in lowering the necessary work W* , 
the nucleation occurs predominantly near the impurity 
or wall. Then the nucleation is called heterogeneous nu- 
cleation. Otherwise it occurs uniformly within the bulk 
and is called homogeneous nucleation. 

Recently, Corti ^ 10 ' 11 and coworkers have studied not 
only the bubble nucleation but the cavity or void forma- 
tion within the stretched liquid using the Monte Carlo 
simulation^, the density functional method*^ and the 
thermodynamic perturbation method— for the Lennard- 
Jones fluid. They conjectured that 

• Conjecture i:The size of the critical cavity is a lower 
bound to the size of critical bubble. 

• Conjecture £:The work of formation of the critical 
cavity is the upper bound to the work of formation 
of the critical bubble. 

where the critical cavity is the largest cavity that can 
be created within the liquid without inducing vaporiza- 
tion^. A larger cavity will immediately drive the sys- 
tem to phase separation into the vapor phase. There- 
fore the critical cavity represents the stability limit of the 
metastable liquid that contains the largest cavity. Based 
on these observation, Corti and coworkers^ have conjec- 
tured that the homogeneous bubble nucleation is induced 
by the homogeneous cavity nucleation. The cavity nucle- 
ation could play an important role to induce subsequent 
bubble nucleation. Then the homogeneous bubble nu- 
cleation could, in fact, be the heterogeneous nucleation 
where the homogenously nucleated cavity would play the 
role of a spherical impurity or surface to induce heteroge- 
neous nucleation of bubbl o 12 ^ 13 ' 14 ! 15 . These conjectures, 
however, are based on the numerical results of the cav- 
ity in the Lennard-Jones system only. Therefore, further 
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confirmation of their findings using a more generic model 
is necessary. 

In this paper, we will use a simple square-gradient 
density functional theory with a triple-parabolic free en- 
ergy^ to study the cavity introduced into the stretched 
liquid 8 911 . The density functional theory (DFT) is 
known to be more reliable than the classical nucle- 
ation theory (CNT) 1,2 in particular near the spin- 
oda l 17 ' 18 ' 19 ' 20 ' 21 . The square gradient approximation 
to DFT is also known to be qualitatively correct not 
only for the liquid with short-ranged interatomic po- 
tentials 2 .^ but also for the liquid with the long-ranged 
potentia l 24 ' 25 . Further approximation using parabolic 
free energy make this square-gradient density functional 
theory more attractive and generic as it does not de- 
pend on the detailed form of interatomic potentials. It 
has elucidated the various aspects of the critical bub- 
ble in stretched liqui d 26 ' 27 ! 28 . For example, we 2 ^ stud- 
ied a scaling rule found for the Lennard-Jones system by 
Shen and Debenedetti 4 - and found that the scaling rule 
is marginally satisfied. Furthermore, this generic model 
can also be used to model the complex liquid such as 
polymer— or even liquid metal**. 

This paper is organized as follows. In Section II we 
present a short review of the triple-parabola model^ 6 - to 
summarize the necessary formula. In Section III, we will 
present the numerical results and discuss the implications 
of the results in light of the conjectures of Punnathanam 
and Corti 9 . Finally Section IV is devoted to the conclu- 
sion. 



II. SQUARE-GRADIENT DENSITY 

FUNCTIONAL MODEL WITH 
TRIPLE- PARABOLIC FREE ENERGY 

A. Triple-Parabolic Free Energy 

Since the model has been explained in detail else- 
where^, we will briefly review the model. In the square- 
gradient density-functional model of the fluid, the free 
energy of the inhomogeneous fluid, such as the critical 
bubble in the stretched liquid is written as 



W 



Aw(0) + c (V0) 2 ) d 3 r, 



(2) 



where = p/pi is the order parameter that represents the 
density p divided by the liquid density /?/. In the triple- 
parabola model of Granasy and Oxtobj*^, the local part 
of the free energy Auj is given by 



Aw(0) 



2 

Ai 

2 



<f>2 — <t>0 



Ap + cj , 



b < (f> A 

!>A < 
(3) 

which consists of three parabola with Ao, A2 > and A x < 
centered at the vapor density 0o, and at the free energy 



barrier 0i (spinodal) , and at the liquid density 02 , which 
we call "vapor", "spinodal" and "liquid" part of the free 
energy. Equation ([2]) correctly represents the necessary 
work of forming bubble or cavity from the uniform fluid 
with relative density 02 since the free energy of uniform 
liquid is given by uj (02) = 0. 

The curvature of parabola Ao and A2 are related to 
the compressibility of vapor and liquid phases, and Ap 
is the free energy difference between the liquid and the 
vapor. Although Ap represents in fact the pressure dif- 
ference AP as Eq. ^ is the grand potential of the open 
system 2 -, we call Ap chemical potential to make the com- 
parison to the previous work 4 *^ easier since AP is pro- 
portional to Ap through AP — p\Ap 2A . We use the 
terminology "over-saturation" when Ap is positive and 
" under-saturation" when Ap is negative. The stretched 
liquid in this study corresponds to the under-saturated 
liquid with negative Ap < 0. 

From the continuity of the free energy Auj(<f) at the 
boundary, the matching densities 4>a and 4>b are given 
by 



Aq0q + I Ai |0! 

A + |Ai| 
A 2 02 + |Ai|0i 

A 2 + |Ai| 



(4) 



while the location 0i of the free energy barrier is given 
by 



(P0O 



h) + ypq 



b 2 )" - 2Ap(p - q) 



, (5) 



with 



A0IA1 



A + |Ai 
A 2 |Ai| 
A 2 + |Ai 



(6) 



The liquid spinodal, where the metastable stretched 
liquid becomes unstable, occurs when the under- 
saturation is given by 



P 



Ap spin = -- (02 - 00 ) 



(7) 



In contrast to the previous models ^ 19 ' 20 ! 21 where the 
compressibility diverges continuously as the spinodal is 
approached (Ap — > A/j, sp i n ), the compressibility remains 
finite until the spinodal point is reached in our triple- 
parabolic model as the curvature Ao and A2 is fixed. In 
the usual 4 model, the divergin g 18 ' 19 ' 29 ' 30 compressibil- 
ity of metastable liquid phase 



< 



A 2 oc 1 



Ap 

Ap spn 



1/2 



(8) 



is predicted. Therefore, the results obtained from this 
triple-parabolic free energy model are somehow artificial 
near the spinodal. 
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FIG. 1: (a) The triple-parabolic free energy from CNT regime 
near the coexistence to the spinodal regime near the liquid 
spinodal for the case-I of section III (Table [XJ . (b) The cor- 
responding critical bubble of homogeneous nucleation at the 
CNT regime (A/i/A/i S pm = 0.3) and at the spinodal regime 
(0.7). 



In Fig. [UJa) we show typical shape of the triple- 
parabolic free energy Atu together with the corresponding 
density profile of the critical bubble in Fig. [Hb). The 
free energy Aw shows a typical double-well form with 
the lower well represents the stable vapor phase and the 
upper well represents the metastable liquid phase. The 
radii r A and tb are the matching radius that satisfies 
4> (, t a) — 4>A and 4> {tb) — 4>b- Since the free energy con- 
sists of three parabolas corresponding to the vapor, spin- 
odal and liquid parts, the density profile of the critical 



bubble also consists of three parts when A/i/A^ s 



0.3 



near the coexistence. However, as the under-saturation 
increases (|A/i| becomes large) and it approaches the liq- 
uid spinodal A/i sp j n , the density profile consists of only 
the two parts corresponding to the spinodal and the liq- 
uid parts as the density never decreases down to the 
vapor density. We use the terminology "CNT regime" 
for the former regime where the classical nucleation the- 
ory (CNT) is expected to be qualitatively correct, and 
"spinodal regime" for the latter where the spinodal nu- 
cleatio n 1 ! 19 ! 20 ! 21 is expected to occur. 



B. Density Profile of Critical Cavity and Bubble 

1. CNT regime 

Density profile of spherically symmetric critical cav- 
ity and bubble can be obtained from the Euler-Lagrange 
(EL) equation SW/5(j> = 0, which corresponds to the sad- 
dle point of the free energy surface. This EL equation 
reads to the differential equation 



d 2 ^ 
dr 2 



r 



0. 



0,1,2 



(9) 



for the three parabolas in Eq. JT]), where Ti — y/\Xi\/2c 
and $i(r) = (f>(r) — fa, and + sign is used for i — 1 and 
— is used for i = 0, 2 for ±. These differential equations 
should be solved with boundary conditions: 



(10) 



together with the condition for the vapor density $o 
given by 

$o (r -> R + ) = -0o (Critical cavity) (11) 
for the critical cavity with the radius R, and 
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dr 



(Critical bubble) 



(12) 



for the critical bubble. The boundary condition for the 
cavity in Eq. (fTTj) implies that the work of formation 
W(R) of the cavity reaches a maximum value as the func- 
tion of the cavity radius R since^L 



— = AttR 2 ^ (R+) k B T = 



(13) 



and 4>(R +) = $ (R + ) + 4>o = 0. This condition can 
also be interpreted as the onset of the drying transition 
of fluid surrounding the cavity^^. The critical cavity 
is not only defined as the saddle point solution of Eq. 
with the largest cavity radius R but corresponds in fact to 
the stability limit as dW/dR > 0. Therefore the critical 
cavity is the limit of the stability of the metastable liq- 
uid that contains the largest cavity. Once the size of the 
cavity exceeds the size of the critical cavity, the stretched 
liquid immediately phase separates into the vapor phase. 
On the other hand, the critical bubble corresponds to 
the intermediate highest energy saddle point of the reac- 
tion of the growth of the nucleus. Once the size of the 
bubble exceeds that of the critical bubble, it increases 
indefinitely and will eventually lead to the phase transi- 
tion of the whole system from the metastable liquid to 
the stable vapor phase. 
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In this CNT regime near the two-phase coexistence, 
both the matching radius and rg exist. When the un- 
der saturation A/j,(< 0) is further increased, the matching 
radius ta approaches zero and disappears for the bubble 
solution (Fig. [I]). Then only the matching radius rs can 
exist. Such a transition never happened to the critical 
cavity as the density always becomes the vapor density 
from the boundary condition Eq. (jlip . In the critical 
cavity, not only the matching radius tb but the radius 
Ta exist up to the spinodal point. 

The solutions of this EL equation Eq. ^ for the crit- 
ical cavity in the liquid (i = 2) and the spinodal (i = 1) 
part of the free energy are given by 

$i(r) = csc(Ti (r A -r B )) 

x {-^ibTb sin (Ti {r~r A )) (14) 
+&1ATA sin (Ti (r - r B )))/r, 

<f> 2 (r) = cf> 2B r B exp (-r 2 r + r 2 r B ) /r, 

while the one in the vapor (i = 0) part is given by 

* (r) = -csch(r (i?-r A )) 

x ($ 0j4 r A sinh(r (r-i?)) (15) 
+0 o i?sinh (r {r~r A )))/r 

The expression for the critical bubble can be obtained by 
setting R = 0. 

Finally, the matching radii ta and tb are determined 
from the simultaneous equation 



(16) 



where $1 are the functions of both ta and Tb- This 
simultaneous equation can be solved numerically us- 
ing standard algorithms such as the Newton-Raphson 
method for both the critical cavity and the bubble. 



2. Spinodal regime 

In this regime near the spinodal, the density of vapor 
bubble remains high even at the center of the bubble. 
Therefore the bubble remains within the liquid well of 
Fig. [T](a), and ta for the critical bubble becomes zero 
(see Fig. [T}. Therefore the i = part of the free energy 
density in Eq. (|3|) and its solution $0 disappears. Then 
the boundary condition Eqs. (fTTj) and (JT3J) for $0 should 
be replaced by the same condition for $1. In fact, such 
a situation can happen only for the homogeneous bubble 
nucleation^ 8 -, but it has never happened to the critical 
cavity as the density always goes down to zero at the 
cavity boundary from Eq. (jlip and the vapor solution 
Eq. (|15p is always necessary to represents cavity density. 
However, for the sake of completeness, we will present 
the cavity solution in the spinodal regime as well. 



d$ 2 




dr 


r=r B ~ dr 


d$i 


d$o 


dr 


, „ dr 



The solution for the EL equation for <i> 2 is the same as 
Eq. I|12p. but the one for $1 now read 



$i(r) = -csc(Ti(i2-rfl)) 

x ($i B r B sin (Tx (r ~ R)) 
+0ii?sin (r x (r-r B ))) 



(17) 



for the critical cavity in the spinodal regime. Again the 
expression for the critical bubble is recovered by setting 
R = 0. 

In this case, the matching radius tb is simply deter- 
mined from the equation 



d$ 2 
dr 



dr 



(18) 



In particular, Eq. (|T8|) for the critical bubble does not de- 
pend on the under-saturation A/i. Therefore, the match- 
ing radius r B is constant in the spinodal regime for the 
critical bubble^. 

In Fig. [IJb), we showed the typical density profiles of 
critical bubble in the CNT regime and in the spinodal 
regime. The critical bubble is larger in the CNT regime 
than that in the spinodal regime. However, the density 
difference between the inside and the outside of the bub- 
ble becomes smaller in the spinodal regime than that in 
the CNT regime. Since the length scale 1/1^ = y / 2c/|A7| 
does not depend on the under-saturation A/i, the shape 
of the critical bubble, in particular, decay length of the 
tail, or the interfacial width look almost constant. 



C. Work of Formation of Critical Cavity and 
Bubble 

1. CNT regime 

Once we know the density profile of bubble and cavity, 
it is straightforward to calculate the work of formation 
W* of the critical cavity and bubble. Since, the results 
for the critical bubble have already been presented else- 
where^, we only show the results for the cavity. To this 
end, we can use the formula 



W* 



4nr z Alu 



1 dAco 



dr 



(19) 



derived by Cahn and HilliarcUi. Using this transforma- 
tion, we can evade the singularity of V</> in Eq. @. Like 
previous author o 8 ' 9 ' 11 , we regard the critical cavity as the 
external hard wall and omit the contribution from r < R 
in Eq. 

Inserting the density profile Eq. (fT2| in the CNT 
regime, we can calculate the integral in Eq. (Til))) ana- 
lytically, and we obtain 



W = W + Wi + W 2 



(20) 
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where 



-in 

Y 



-Afi- 



2 - <P0 
(-$ 1A r A + QlBTB 
„2 



A/t + UlQ 



w 2 



2ttJAi^ 

r? 

- ri ($i A r 2 A + $i B r 2 B ) cot (r t (r A - r B )) 

- Ti (® 1A + r A r B csc (Tx (r A - r B ))) 
2n\2(j)2®2Br B 



(21) 



r 2 

1 2 



(i + r B r 2 



for both the critical cavity and critical bubble, while 
VT = — r\Afi + ° v (7T i? + Qoata 

1 

+r (0 o i? 2 + $ 0A r^) coth (To (R - r A )) (22) 
+r (0o - $oa) M^csch (r [R - r A ))) 

for the critical cavity. The expression of Wo for the crit- 
ical bubble is recovered when R = 0. Note that A/i < 
for the critical bubble in the under-saturated liquid. 



2. Spinodal regime 

Since, we only have the solution <J> 2 in Eq. (TTJ} and $i 
given by Eq. (fT7|) for the critical bubble in the spinodal 
regime, we have 

W = W x + W 2 , (23) 
where W 2 is given by Eq. (f2"Tj) but W\ is given by 

47T 



2vr|Ai|0i 



A/i + LUq 



(24) 



+r x (-0ii? 2 + $i B r|) cot (r x (JE - r B )) 
+ri (0i - $ib) r B i?csc (Ti (iE - r fl ))) 

for the critical cavity. The expression of W\ for the criti- 
cal bubble is recovered when R = 0. However, this equa- 
tion Eq. (|2"3")l will not be used for the critical cavity as it 
corresponds to the unphysical solution with higher free 
energy. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In order to check the validity of the two conjectures 
by Punnathanam and Corti 9 mentioned in the introduc- 
tion derived from the numerical results from the Lennard- 
Jones fluid, we use this generic square-gradient density- 
functional model with triple-parabolic free energy and 
compare the various properties of the critical cavity and 
bubble. We use several typical free energy parameters 



TABLE I: Three sets of free energy parameters used in this 
work to study the various properties of critical cavity and 
bubble. 



model 


c 


00 


02 


A 


Ai 


A 2 


case-I 


1.0 


0.1 


1.0 


1.0 


-0.5 


0.3 


case-II 


1.0 


0.1 


1.0 


1.0 


-1.0 


0.7 


case-Ill 


1.0 


0.1 


1.0 


1.0 


-2.0 


0.9 



used before 2 ^ to confirm the results further. The three 
sets of the free energy parameters used are summarized 
in Tabled] 

Figure [2] compares the matching radius r A and r B as 
well as the critical radius R of the critical cavity and 
bubble^ as the functions of the scaled under-saturation 
A/t/A/t S p for the case-I to III. The critical radius is the 
maximum radius R for which the saddle point solution 
for Eq. © can exist. We manually increased the radius 
R of the cavity and monitor the stability of the solution 
and determined the maximum radius R. 

Figure [2] also shows the matching radius r A and r B 
of the critical bubble reported before 2 ^ as the functions 
of the scaled under-saturation A/t/A/i sp together with 
the radius R of the critical cavity for the case-I to III. 
The matching radius r A becomes zero and the density 
profile is represented by Eq. (flT)) as the under-saturation 
enters the spinodal regime, while the matching radius r B 
survives and it becomes constant and independent of the 
under-saturation in the spinodal regime as predicted from 
Eq. (|17p . We also show the equimolar dividing radius r B 
of the critical bubble determined from 28 



(0(r) — 0or g ) 47rr dr = (0 2 — 0(r)) 47rr dr, 

10 Jm 

(25) 

where 



0(r -> 0) 



(26) 



is the density at the origin of the bubble. The divid- 
ing radius r B is constant as the matching radius r B in 
the spinodal regime. In fact, these constant radii are 
due to the fact that the parameter A 2 which represents 
the compressibility of the metastable liquid remains con- 
stant. Therefore the constant radius in the spinodal 
regime up to the spinodal is merely the artifact of our 
model. Rather, the diverging radius is expected as the 
compressibility of metastable liquid diverges at the spin- 
odaUi d 9 ' 20 according to Eq. ([8]). In contrast, the critical 
radius R of the critical cavity converges to zero at the 
spinodal. 

In contrast to the critical bubble, the matching radius 
r A of the critical cavity does not disappear in the spin- 
odal region but decrease monotonically and approaches 
zero at the spinodal. This is due to the boundary con- 
dition Eq. (jTTJ) and the fact that the density always de- 
creases down to a vacuum. The density profile of critical 
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FIG. 2: The matching radii ta (circle) and tb (triangle) as 
well as the critical radius R (square) of the critical cavity. 
Also shown are the matching radii ta (solid curve) and tb 
(chain curve) as well as the dividing radius ro (dashed curve) 
of the critical bubble for the case-I to III. 



cavity remains within the CNT regime down to the spin- 
odal. Since the matching radius ta approaches zero at 
the spinodal, the critical radius R also approaches zero 
as R should always be smaller than the matching ra- 
dius ta of the critical cavity. The matching radius Tb 
remains finite even at the spinodal, which indicates the 
finite width of the liquid-vapor interface of the cavity. 
No matter what the definition of the size of the critical 
cavity and bubble, the size of critical cavity represented 
by R or rs seems the lower bound of the corresponding 
radius such as the dividing radius rjj and the matching 
radius re of the critical bubble as shown in Fig. O 

In fact, Figure [3] clearly indicates that not only the 
radius R but the overall profile of the critical cavity is 
smaller than the critical bubble. In the CNT regime of 
the low under-saturation, in particular, the density pro- 
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FIG. 3: (a) The density profile (order parameter) of the 
critical cavity compared with that of the critical bubble in the 
CNT regime of the low under-saturation A^/A^ sp i n = 0.15 
for the Case-I. (b) The density profiles in the spinodal regime 
of the high under-saturation A/i/A/i sp i n = 0.90. 



file of critical cavity is similar to the corresponding crit- 
ical bubble. Not only the size but also the width of the 
liquid-vapor interface is similar magnitude. However, the 
liquid within the critical bubble is more pushed outward 
compared to the liquid within the critical cavity even 
though the density becomes zero for r < R for the crit- 
ical cavity (Fig. Ela)). The typical radius of the bubble 
ro or rg is always larger than R or tb for the critical 
cavity as shown in Fig. [2j 

The situation is less clear in the spinodal regime as 
the density profile of the critical cavity and bubble are 
totally different (Fig. 02b)). The density of cavity goes 
to zero at the critical radius R while the density of the 
critical bubble within the bubble is almost the same as 
that of the liquid. Therefore the width of the liquid- vapor 
interface is also different, and the interface of bubble is 
much wider than that of the cavity. Yet, we can say that 
the cavity is smaller than the bubble as the typical size 
of the cavity, for example, is always smaller than that 
of the bubble. In particular, not only the typical size ta 
but the critical radius R for the critical cavity approaches 
zero at the spinodal, while the dividing radius rjj for the 
critical bubble remains finite^. 

The density profile of critical cavity is the saddle point 
solution of the EL equation Eq. ([9|) with the lowest free 
energy and maximum radius R. As we increase the radius 
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of cavity, the work of formation W increases monoton- 
ically because dW/dR > from Eq. (fT3|) . As soon as 
the radius of the cavity exceeds the critical radius R that 
corresponds to dW/dR = 0, the free energy jumps to the 
higher energy "excited-state" of Eq. ([9]) and the oscil- 
lating density profile appears. The origin of this oscilla- 
tion is not the same as the one observed in the previous 
Monte Carlo simulation^ as the latter comes from the 
excluded volume effect. Our square-gradient approxima- 
tion cannot include such an excluded volume effect, and 
our oscillatory solutions are not real. Rather such an os- 
cillatory solution indicates that the stable solution which 
smoothly connects to the liquid density at infinity cannot 
exisl£. The cavity with a radius larger than the critical 
radius R induces more density oscillation and a higher 
free energy. Therefore, we interpreted the appearance of 
oscillatory solution as the limit of stability of the cavity. 
The solution with maximum radius R without density 
oscillation corresponds to the critical cavity which is the 
cavity of the maximum size that can be accommodated 
into the stretched fluid. 

The work of formation W* of the critical bubble and 
cavity are compared in Fig.[4j The work of formation for 
the cavity is always much larger than the critical bub- 
ble. The work of formation of critical cavity approaches 
zero at the spinodal. In contrast, the work of formation 
of critical bubble calculated from the CNT erroneously 
remains finit e 23 ' 29 even at the spinodal. Although both 
the work of formation of cavity and bubble from the DFT 
approaches zero, the former seems always larger than the 
latter in accordance with the results from a more sophis- 
ticated non-local DFT for the Lennard- Jones fluid 9 . 

From the results for the case-I to III in Figs. [5] and 
IU it seems fair to say that the size of the critical cav- 
ity is always smaller than the size of the critical bubble, 
while the work of formation of the critical cavity is always 
larger than that of the critical bubble in accordance with 
the conjecture 1 and conjecture 2 of Punnathanam and 
Corti 9 - mentioned in the introduction. In particular, the 
work of formation for the critical cavity is much higher 
than that of the critical bubble. 

Figure [5] shows the work of formation of the critical 
cavity and the bubble as the function of the cavity radius 
R of the critical cavity or the dividing radius t\d of the 
critical bubble. The work of formation W* for the critical 
cavity scales with the parameter R, while the one for the 
critical bubble does not scale with the dividing radius td- 

In the previous paper—, we have checked the scaling 
rule for the various quantities of the critical bubble of 
our square-gradient DFT-model as a function of scaled 
under saturation Afi/ A/i sp ; n . In contrast to the Lennard- 
Jones fluid where such a scaling rule is indeed valid^, we 
found that although some quantities show almost perfect 
scaling relations near the spinodal, the work of forma- 
tion divided by the value deduced from the CNT shows 
no scaling 28 . Similar scaling rules for cavity are also 
found in the Lennard- Jones fluid by Punnathanam and 
Corti 9 -. In order to check the scaling rule of cavity in our 
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FIG. 4: The work of formation W* of the critical cavity 
(circle) compared with that of the critical bubble calculated 
from DFT (solid line) as well as that from CNT (dashed line) 
for the case-I to III. 



DFT-model we have plotted the critical radius R and the 
matching radius rs of the critical cavity divided by rs 
of the critical bubble as the function of scaled undersat- 
uration A/x/A/i sp i n in Fig. [HIa). The critical radius R 
as well as the matching radius rs of cavity show sim- 
ilar scaling relations in particular near the spinodal to 
the various radii of critical bubble 28 of the same DFT- 
model. However, Fig.[Sfb) clearly indicates that the work 
of formation W* of the critical cavity as the function of 
the scaled undersaturation A/i/A/x sp i n shows no scaling 
again2£. 

We should note that the cavity and bubble in Fig. 
represents critical cavity and bubble exactly at the saddle 
point as we consider the open system within the frame- 
work of the grand canonical ensemble. Therefore this 
diagram is not for the work of formation W(n, r) of cav- 
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FIG. 5: The work of formation W* of the critical cavity 
(circle) and bubble (solid line) as the function of the cavity 
radius R or the dividing radius ro of the bubble. 
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ity or bubble with arbitrary number of molecule n with 
arbitrary size r or R in the canonical ensemble2*ii. Al- 
though this diagram is for the critical cavity or bubble, 
it is expected that the work of formation of any cavity 
or bubble with arbitrary size in the canonical ensemble 
would show a similar trend as the function of its radius 
R. Then, it is expected that the (sub-critical) cavity, 
which is smaller than the critical cavity and is not large 
enough to induce instability, would always exist near the 
spinodal as both the work of formation W* and the size 
of the critical cavity R of even the largest critical cavity 
approaches zero at the spinodal. 

In our DFT-model, the size of the critical bubble ap- 
proaches finite minimum size while the work of formation 
approaches zero at the spinodal—. On the other hand, 
both the size and the work of formation of the critical 
cavity approaches zero at the spinodal. Even though the 
work of formation of critical cavity is higher than that of 
the critical bubble as shown in Figs. |3J the sub-critical 
cavities smaller than the critical cavity could be easier to 
form than the much larger and diffuse critical bubble near 
the spinodal. Then the bubble formation might be pre- 
ceded by the (sub-critical) cavity formation at least near 
the spinodal, and the homogeneous nucleation of bubble 
might be in fact the heterogeneous nucleation where the 
homogeneous bubble nucleation starts from the cavity 
formation. As the size of the sub-critical cavity grows, it 
transforms into the bubble by accommodating the liquid 
molecule into the cavity and lowers its free energy (work 
of formation) before reaching the critical cavity. In such 
a case the phase separation of stretched liquid to vapor 
occurs via the gradual heterogeneous bubble nucleation 
around the sub-critical cavity rather than the explosive 
critical cavity formation that corresponds to the stability 
limit of metastablc liquid. In fact, recent molecular dy- 
namics 32 and Monte Carlo! studies of bubble formation 
and boiling show that the bubble nucleation is initiated 
by the spontaneous formation of vacuum cavity. 

In the previous paper, we have pointed out that 
this triple-parabolic free energy has an artifact that the 
isothermal compressibility remains constant as the pa- 
rameter A, is fixed in Eq. @, which reflects in the finite 
and constant size of the critical bubble represented by the 
constant dividing radius rjj or the matching radius Tb- 
The same behavior is expected when we use the double- 
parabola model 2 ^. It is expected, in fact, the isothermal 
compressibility k of the liquid phase diverges as^ 



K oc A, oc 1 



A^ S pi r 



-1/2 



(27) 



from Eq. ([5]) of the quartic (f) A field models. Then the 
size of the critical bubble should diverge a o 17 ' 18 i 29 



FIG. 6: (a)Scaled radius R and the dividing radius tb of 
critical cavity as the function of the scaled undersaturation 
A/i/A/i S pi n . (b) Scaled work of formation W* of critical cavity 
as the function of the scaled undersaturation A^i/A/i sp m. 
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-1/4 



(28) 



when A/i/A/i sp j n — > 1 as the density profile of the crit- 
ical bubble is given by Eq. (jTTJ) with R = and rs is 
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determined from in the spinodal regime. Then, the 
story for the critical bubble in real materials would be 
slightly different from that of our triple-parabolic model, 
in particular, in the spinodal regime. On the other hand, 
since the density profile of the critical cavity is always 
given by Eqs. (fl4|) and (fl~5|) of the CNT regime, the size 
of the critical cavity R will be relatively insensitive to the 
A 2 . 

It is clear from Fig. [2] and Eq. (|28| that there exists a 
minimum size of the critical bubble, while the size of the 
critical cavity approaches zero. Furthermore, the critical 
size R of the cavity is smaller than the size of the crit- 
ical bubble represented by the dividing radius Tb from 
Fig. E^a) . Even though the critical size of the cavity is 
smaller than the bubble, it is also clear from Fig. E^b) 
that the work of formation of the critical cavity is always 
larger than that of the critical bubble. 

Since both the size of the critical cavity R and the work 
of formation W* approaches zero as we move toward the 
spinodal, it is expected that the sub-critical cavity rather 
than much larger critical bubble could easily be formed 
and induce phase transformation from metastable liquid 
to vapor near the spinodal. The previous picture of ho- 
mogeneous nucleation assumes that the long-wavelength 
fluctuation induces the critical bubble formation. Near 
the spinodal, such a homogeneous nucleation called spin- 
odal nucleation would occur cooperatively and the sud- 
den and the rapid phase transformation could occur as 
the system becomes nearly unstable. Our simple calcu- 
lation using DFT model indicates that this scenario of 
phase transition near the (mean-field) spinodal may not 
be true and the sub-critical cavity rather than the critical 
bubble could play some role to induce the rapid spinodal 
nucleation. 



IV. CONCLUSION 

In this work, we have studied the size and the work 
of formation of the critical cavity using a generic square- 
gradient density-functional model with triple-parabolic 
free energy. We have used this model as it is generic 
in the sense that it does not depend on the interatomic 
potential and it has already been used to study the scal- 
ing properties of the critical bubble— . We pay partic- 
ular attention to the comparison of the critical cavity 



with the critical bubble of the homogeneous nucleation 
to check the conjecture made by Corti and coworker— for 
the Lennard- Jones fluid. 

Aside from the fact that our DFT model assumes con- 
stant compressibilities A,, i = 0, 1, 2 in Eq. ([3]) and a con- 
stant square-gradient coefficient c in Eq. it seems 
certain from our numerical results that the two conjec- 
ture made by Punnathanam and Corti 9 - on the size and 
the work of formation of the critical cavity are valid. 
From the behavior of the critical cavity near the spin- 
odal, sub-critical small cavity seems easier to form than 
the critical bubble in stretched fluid with negative pres- 
sure near the spinodal. This conclusion for the sub- 
critical cavity and critical bubble near the spinodal re- 
main correct even if we include the diverging isother- 
mal compressibility at the spinodal predicted from the 
4 -field theory ^ 19 ' 29 . Then it could be argued that the 
sub-critical cavity 9 -^ plays some role to induce critical 
bubble nucleation in the stretched fluid at least near the 
spinodal. If the sub-critical cavity plays the crucial role 
in liquid-to- vapor phase transformation (vaporization), 
the real picture of the phase transformation should be 
more complex. Then the real picture of phase transfor- 
mation is totally different from the traditional picture 
based on the homogeneous nucleation where the phase 
transformation occurs through the formation of the crit- 
ical bubble from the fluctuation and subsequent growth 
of the bubble. It would be, however, difficult experimen- 
tally to judge whether the phase transformation near the 
spinodal proceeds via the cavity or bubble as the liquid 
phase becomes nearly unstable near the spinodal where 
the fluctuation plays the dominant role and smears the 
detail of the phase transformation. 

Finally, the dynamics of cavitations would be more in- 
teresting for practical purposes. Various numerical meth- 
ods such as the Monte Carlo^ 3 ., the molecular dynam- 
ics^ and the lattice Boltzmann method 3 ^ have been de- 
veloped to study the dynamic of cavity or bubble for- 
mation and the evaporation. Recently, we have devel- 
oped a numerical method based on the time dependent 
Ginzburg-Landau model combined with the cell dynam- 
ics method to study the dynamics of nucleation in various 
situation o 36 ! 37 . It will be interesting to use this cell dy- 
namics method to study the dynamics of cavity or bubble 
formation. 
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